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Two classes of models of driven disordered systems that exhibit history-dependent dynam- 
ics are discussed. The first class incorporates local inertia in the dynamics via nonmonotonic 
stress transfer between adjacent degrees of freedom. The second class allows for proliferation 
of topological defects due to the interplay of strong disorder and drive. In mean field theory 
both models exhibit a tricritical point as a function of disorder strength. At weak disorder 
depinning is continuous and the sliding state is unique. At strong disorder depinning is 
discontinuous and hysteretic. 

I. INTRODUCTION 

Nonequilibrium depinning transitions from static to moving states underlie the physics of a 
wide range of phenomena, from fracture propagation in heterogeneous solids to flux flow in type-II 
superconductors One class of models, overdamped elastic media pulled by a uniform force, has 
been studied extensively. These exhibit a nonequilibrium phase transition from a pinned to a sliding 
state at a critical value Ft of the driving force. This nonequilibrium transition displays universal 
critical behavior as in equilibrium continuous transitions, with the medium's mean velocity v 
acting as an order parameter P, Q| . In overdamped elastic media, the sliding state is unique and 
no hysteresis can occur Q|. 

The elastic medium model is often inadequate to describe real physical systems that exhibit 
plastic response on various scales. Plasticity, used here in a loose sense, may arise in systems with 
an underlying periodic structure, such as vortex lattices or charge density waves, for strong disorder. 
This yields large deformations of the driven medium with the proliferation of topological defects, 
which are continuously generated and healed by the interplay of drive, disorder and interactions. ^, 
In other situations, such as crack propagation in heterogeneous solids P, and the motion of a 



helium-4 droplet contact line on a rough surface |9|], the dissipative elastic medium model is made 
inadequate by the presence of inertial effects or other nonlocal stress propagation mechanisms. In 
both cases the depinning transition may become discontinuous and sometime hysteretic. 

Several coarse-grained models of driven extended systems that can lead to history-dependent dy- 



namics have been proposed in the literature 



E, 10, U, 12, 3 



Here we examine two classes of such 
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models. In the first class of models the displacement of the driven medium from some undeformed 
reference configuration remains single- valued, as appropriate for systems without topological de- 
fects, but the elastic interactions are modified by assuming local underdamped dynamics. This 
is modeled via a linear stress-strain relation, where the stress transfer between displacements of 
different parts of the manifold is nonmonotonic in time. Models of this type have been used to 
incorporate the effect of inertia or elastic stress overshoot in crack propagation in solids [l^l . 14]. 
A related model was proposed by the author and collaborators as an effective description of topo- 
logical defects in the driven system 
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15. 



lfj. Below we show that the presence of topological 
defects in a solid can be described as a viscous force that allows a moving portion of the medium to 
overshoot a static configuration before relaxing back to it, so that these two models are identical. 



In the second class of models topologica 
that displacements are single- valued [l 



defects are explicitly allowed by removing the constrain 



12j. We consider a simple realization of this type of mod- 



els proposed some time ago in the context of charge density waves, where the scalar displacement 
describes the phase of the electronic condensate and the system only exhibits one-dimensional pe- 
riodicity along the direction of motion. It is obtained by assuming a nonlinear coupling among 
neighboring displacements that incorporates the crucial feature that the displacement becomes un- 
defined at the location where the amplitude of the order parameter collapses. This is incorporated 
in the model as an instantaneous jump of the displacement or phase of an amount equal to its 
period, knows as phase slip. 

The two classes of models exhibit remarkably similar behavior in the mean-field limit, where 
the interaction are assumed to have infinite range. As a function of disorder strength, they both 
have a tricritical point separating continuous from discontinuous depinning transitions. Depinning 
is continuous at weak disorder with mean-field critical exponents that are identical to those ob- 
tained for a dissipative elastic medium. Above a critical disordered strength, depinning becomes 
discontinuous and hysteretic. The origin of the hysteresis is, however, different in the two models, 
as discussed below. Preliminary numerical work also suggests that differences may arise in finite 
dimensions. 



II. TWO CLASSES OF MODELS: VISCOUS AND PHASE SLIP COUPLINGS 

We restrict ourselves to the dynamics of a scalar field describing deformations along the direction 
of mean motion and discretize space, denoting by Ui(t) the displacement of the i- degree of freedom 
from some undeformed reference configuration, at time t. We stress that all models discussed 
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are coarse-grained ones, with u, representing the displacement of a region pinned collectively by 
disorder. The microscopic dynamics is always assumed to be overdamped, but velocity-dependent 
couplings can arise upon coarse-graining. The equation of motion in the laboratory frame is written 
by balancing all the forces acting on each segment of the driven medium as 

d t u i = a a ({u i },t)+F + f p (i,u i ) , (1) 

where time has been scaled so that the friction coefficient is unity, F is the external driving force, 
f p is the pinning force, and a a represents the stresses due to interactions with the neighbors. The 
label a will be used below to identify two ways of modeling the coupling. For driven periodic 
manifolds, such as vortex lattices or charge density waves, the pinning force is periodic in the 
displacement Uj, with < ui < 1, and can be written as f p (i,Ui) = hiYiui — 7j), where Y(u) is a 
periodic function specified below, hi are random forces with distribution p{hi), and 7, are random 
phases chosen independently and uniformly in [0, 1). 

First, we consider a class of models with single- valued displacements Ui and a linear stress-strain 
relation of the general form, 

ri 



a ({ Ui },t) = [ dt' , (2) 

where the sum is over the Z sites j that are nearest neighbors of i. A suitable choice of the stress- 
transfer or memory function, Jj"(i), yields various models described in the literature. The familiar 
elastic model is obtained by assuming instantaneous stress transfer, i.e., JE(t) = K5(t), with K an 
elastic constant. More generically, all monotonic models, defined as those with Jij(t) > 0, for all 
i,j, t, exhibit a continuous depinning transition with universal critical behavior and a unique sliding 
state The choice jY = r]dt5(t) yields purely viscous stresses, i.e., <t v = v12{j) ~~ u i(t)]> 

provided we identify Vi = iii = dtUi with the local flow velocity of the medium. In this case 
depinning always occurs at F = for pinning force distributions without a finite lower bound, 
but there is a critical point above which the system can switch discontinuously and hysteretically 
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la ] . The author and collaborators 



between a macroscopic slow moving and a fast moving state 
recently considered models where the stress transfer function is taken to have the form appropriate 
for a viscoelastic fluid, that responds elastically on short time scales and flows viscously at long 
times This viscoelastic coupling was proposed as an effective way of describing the local slip 
due to dislocations generated at the boundaries between coarse-grained degrees of freedom. The 
connection between the viscoelastic model and the presence of free dislocations in the medium was 
made more precise in Ref. where it was shown that the equations describing the dynamics of 
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equilibrium deformations of a two-dimensional lattice with a finite concentration nd of f ree annealed 
dislocations can be recast in the form of the phenomenological equations of a viscoelastic fluid 
introduced many years ago by Maxwell j^jj. In a scalar version of the equations of viscoelasticity, 
local compressional stresses are written in the form given in Eq. (J2J), with 

J™{t) = Koo S(t) - K °°~ K \ -t/T (3) 

r 

where and Kq are the high and low frequency compressional moduli, respectively, and r ~ 1 /rid 
is a microscopic relaxation time j^fll ]. Shear stresses have a similar form, with shear moduli replacing 
compressional ones and Go = as a fluid has no zero-frequency elastic restoring forces in response 
to shear stresses. The nonzero long-wavelength compressional elasticity {Kq ^ 0) is associated 
with density conservation and plays a crucial role in controlling the physics of depinning in driven 
lattices. For this reason in our scalar model we assume a coupling of the form © in all directions. 
On time scales t short compared to t, the contribution to the stress coming from the second term 
in Eq. Q is negligible compared to the first and <7 VE reduces to the stress of an elastic solid, 

a™{{ Ui },t « r) « ^£M*) "«<(*)] • (4) 

(?) 

Conversely, for t » r, one can expand the relative displacements for t' ~ t, and 

a VE {{ Ui }, t » r ) « - 5>j(i) - Ui{t)] + | 5>i(*) - rki(i)] , (5) 

<?> <?> 

with i'T = Kq and ?7 = {K^ — Kq)t. The model of driven depinning studied below is obtained from 
Eq. Q with the simplified form given in Eq. © for the stress strain relation, 

dm = ^ J>j(t) - «i(t)] + | ^fe(t) - + F + / p (i,ni) , (6) 
(?) (?) 

and will be referred to as viscous/elastic model (VE). The presence dislocations generated by 
disorder is incorporated in a mean-field-type approximation as a local inertial response of the 
driven system embodied by the viscous coupling of strength r/ ~ l/n^. This model assumes a fixed 
number of topological defects and does not describe the creation and annihilation of dislocations 
due to the interplay of drive, disorder and interactions. Furthermore, in a driven disordered solid 
unbound dislocations can be pinned by disorder and do not equilibrate with the lattice. The 
resulting dynamics cannot be described as a locally underdamped response of the systems. The 
model given in Eq. (jHJ) may be used to describe some of the effect of topological defects near 
depinning, but becomes inapplicable at large driving forces where dislocations recombine as the 
lattice reorders. 
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14j |. although the random pinning 



Before studying the dynamics of the VE model, it is useful to discuss its relationship to other 
models studied in the literature. In particular, the form of the stress transfer function given in 
Eq. © is the same as that used recently by Fisher and Schwarz to incorporate the effect of stress 
overshoot on propagation of cracks in heterogeneous solids 
force considered there is not periodic. These authors consider an automaton model where time is 
discrete. It is straightforward, to define an automaton version of our VE model, where both the 
displacement Ui and time are discrete. The displacement takes integer values and the automaton 
is updated according to the rule 

Ui {t + 1) = Ui + Q{Fi( Ui )) , 

v i (t + l) = @(F i (u i )), (7) 
where Vi can have values and 1 and Fi{u{) is the total force at site i given by 

Fi(ui) = ^J2 M*) ~ "«(*)] + Z (l + v) S v i® + F+ ~ fp{€) ■ (8) 



The pinning force f p (i) becomes a random number chosen uniformly from an interval [0, ho] 
The automaton can be obtained in the limit of very deep periodic pinning wells, when the dynamics 
is dominated by discrete events corresponding to jumps of the displacement from one well to 
the next. As long as v « 1, where 1 is the artificial upper limit of the mean velocity in the 
discrete model, the automaton dynamics mimics the continuous time dynamics reasonably well. 
The automaton version of the viscoeleastic model given in Eqs. ((7J) and © is identical its dynamics 
to the model of crack propagation with stress overshoot studied by Schwarz and Fisher, provided 
the strength M of the stress overshoot is identified with the combination 77/(1+77). The two models 
differ in the type of pinning considered as the random force f p (i,Ui) used in Refs. 13, ^| is not 
periodic. By establishing the connection between these two models we have shown that distinct 
physical mechanisms (inertia, nonlocal stress propagation, unbound topological defects) at play in 
different physical systems can be described generically by a coarse-grained model that includes a 
coupling to local velocities of the driven manifold. 

In the second class of models topological defects are explicitly allowed by removing the constraint 
of single-valued displacements. At a strong pinning center, deformations of the driven medium can 
be large and lead to the accumulation of a large strain. When the distortion is released through 
a collapse of the amplitude of the order parameter (i.e., the creation of a topological defect), the 
displacement abruptly advances of an amount of order 1, while the amplitude quickly regenerates. 
This process is known as phase slippage in superconductors and superfluids. On time scales large 



compared to those of the microscopic dynamics, it can be described approximately as a "phase slip" : 
an instantaneous (modulo 1) hop of the displacement of an integer unit, modeled as a coupling 
periodic in the difference in displacements between neighboring degrees of freedom, 

a PS ({u t },t) = J J>n27r[ % (t) - Ui (t)} . (9) 

ti) 

If the relative displacements are small, this reduces to an elastic coupling of strength K. This model 
has been studied before in the mean-field limit by Strogatz and collaborators for a sinusoidal pinning 
force In this case depinning is always hysteretic in mean field. Very recently we were able to 
solve the model in mean field for arbitrary pinning potential and show that more generically both 
continuous and discontinuous depinning transitions are obtained as the parameter K is varied • 
Below we compare the behavior of these two classes of models: models where the displacement 
remains single valued and deviations from elasticity are introduced via local inertial couplings and 
models where the displacement ceases to be single valued and topological defects can be generated 
for strong disorder. 



III. MEAN-FIELD SOLUTION 



The mean-field approximation for the VE model is obtained in the limit of infinite-range elastic 
and viscous interactions. Each displacement then couples to all others only through the mean 
velocity, v = and the mean displacement, u = iV" 1 We look for solutions 

moving with stationary velocity, so that u = vt. Since all displacements U{ are coupled, they 
can now be indexed by their disorder parameters 7 and h, rather than the spatial index i. The 
mean-field dynamics is governed by the equation 

(1 + rj)u{h, 7) = K(vt - u) + F + rjv + f p (u; h, 7). (10) 

It is useful to first review the case i] = 0, where Eq. (jlUj) reduces to the mean field theory of a driven 
elastic medium worked out by Fisher and collaborators Q]. The mean field velocity is determined 
by the self-consistency condition (u(h, 7) — vt)h^ = 0, where the subscripts h,j indicate averages 
over the distribution of pinning strengths, p(h), and over the uniformly-distributed phases, 7. 
For piecewise harmonic pinning, Y(u) = 1/2 — u, for < u < 1, no moving solution exists for 
F < Ft = ( 2(K+h) } h ' J us t & bove threshold the mean velocity has a universal dependence on the 
driving force, with v ~ (F — Fxy- In mean-field the critical exponent (3 depends on the shape of 
the pinning force: (3 = 1 for the discontinuous piecewise harmonic force and /? = 3/2 for generic 



smooth forces. Using a functional RG expansion in 4 — e dimensions, Narayan and Fisher 
showed that the discontinuous force captures a crucial intrinsic discontinuity of the large scale, 
low-frequency dynamics, giving the gener al result = 1- e/6 + C(e 2 ), in reasonable agreement 
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26(| . For simplicity and to reflect the 



with numerical simulations in two and three dimensions 
"jerkiness" of the motion in finite-dimensional systems at low velocities, we use piecewise harmonic 
pinning. 
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FIG. 1: Mean-field solution of the VE model with a piecewise parabolic pinning potential, p(h) = 8(h — ho) 
and T) — 5. Left frame: phase diagram. Here "coexistence" refers to multistability of the solutions to 
the equations of motion. Right frame: velocity versus drive for ho/K = 0.5 (blue), ho/K = 1 (red) and 
ho/K = 2 (black). Also shown for ho/K = 2 are the discontinuous hysteretic jumps of the velocity obtained 
when F is ramped up and down adiabatically. 



When 7] > 0, the nature of the depinning differs qualitatively from the r] = case in that 
both unique and multi-valued solutions can exist depending on the values of the parameters. The 
solution for general rj can be found from that for rj = by substituting the effective driving force 
G = F + r]v for F in the v^ = o(F) relation and scaling the velocity down by 1 + r/. The linear 
transformation F = G — rjv then gives the general v (F) curve. The mean velocity in the sliding 
state is given by the solution of 

F — Ft = v[l — M( V , K)\ + < — - h)[ J +h) f {l+v)v _ ^ h , (H) 

where M(rj, K) = (1 + r])( ij^^v^ h- For M(rj, K) < 1 there is a unique sliding solution with mean 
velocity near threshold given by 
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with (3 = 1. The condition M(j],K) = 1 determines a critical line separating unique from multi- 
valued solutions v(F). The phase diagram in the (F, ho) plane is shown in Fig.nfor p{h) = 5{h—ho) 
(provided if ^0, the topology of the phase diagram does not depend on the form of p(h)). There 



is a tricritical point at (h c ,F c = Ft), with h c = K/(\/l + W— 1)- For ho < h c , a continuous 
depinning transition at Ft separates a stable pinned state [2j| from a sliding state with unique 
velocity given by Eq. (|12jl. In finite dimensions, we expect this transition to remain in the same 
universality class as the depinning of an elastic medium (77 = 0). This is corroborated by numerical 
studies and analysis by Schwarz and Fisher of the related stress overshoot model (but with 
non-periodic pinning). For ho > h c The v(F) curves are multivalued, which leads to hysteresis 
when F is ramped up and down adiabatically. The hysteresis is easily understood as a consequence 
of the global nature of the viscous coupling in mean field, where the driving force is replaced by 
an effective drive F + rjv. Clearly this does not affect the static state where v = 0, so that upon 
ramping up the driving force from the pinned state the system always depins at Ft = Ft- When 
the force is ramped down from the sliding state where v 7^ 0, the system sees a larger effective 
drive F + rjv and repins at the lower value It is important to appreciate the crucial role of 
a finite value of K in Eq. (JBJ). When K = each degree of freedom has its own velocity Vi and 
for any broad p(h) there are always some de gree s of freedom that experience zero pinning force, 



yielding no stable pinned phase at any F > 



12] . For finite long-time elasticity, i.e., when K 7^ 0, 



the elastic forces or particle conservation enforce a uniform time-averaged velocity for all degrees 
of freedom and one obtains a stable pinned phase for F < Ft- Finally, we note that the VE model 



is also closely related to a model of sliding CDWs that 
normal carriers by adding a global velocity coupling [11 



incorporates the coupling of the CDW to 
2^] to the Fukuyama-Lee-Rice model. 



We now turn to the phase slip model, where the stress is given by Eq. @. In mean field 
theory the nonequilibrium state can be described in terms of two order parameters. The first is 
the coherence of the phases, measured by the amplitude r of a complex order parameter defined 
by 

N 



i=l 



^X> i2 ™s (13) 



with tp a mean phase. In the absence of interactions among the phases or external drive, the Ui 
are locked to the random phases 7, and the state is incoherent, with r = 0. In the opposite limit of 
very strong interactions we expect perfect coherence of the static state, with all Ui becoming equal 
to the mean phase and r — > 1. The second order parameter is the average velocity of the system, 
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given by 



1 N 



(14) 



i=l 



The mean velocity is the order parameter for the transition between static and moving phases. In 
the sliding state ip = vt. The mean-field equation of motion is 



u(h, 7) = F — Kr sin(27ru — ip) + hY(u — 7) 



(15) 



where the effective coupling is proportional to the coherence r. The self-consistency condition for 
the mean field theory is given by 



re 



id' 



dj / dhp(h)e i2nu( - h ^ 



(16) 



In we developed a general method for solving the mean field equations (|15|) and (|16ft for 
arbitrary pinning potential. The resulting phase diagram for a piecewise parabolic pinning potential 
is shown in Fig. El We see immediately that this is richer than the phase diagram for the VE model, 
as now both disorder-driven transitions between static coherent and incoherent phases as well as 

n 

force-driven transitions between static and moving phases can occur. Although, as shown in jl 71 ] . 
the detailed shape of the phase boundaries depends on the specific form of the pinning potential in 
mean field, the types of phases and the schematic topology of the phase diagram are general. At 
low driving forces both coherent (CS) and incoherent (IS) static phases are stable. The transition 




FIG. 2: Mean-field phase diagram for the PS model with a piecewise parabolic pinning potential and 
p(h) = S(h — ho). In the cyan region both stable CM and IS solutions exist. In the small purple region 
near the tricritical point hysteretic depinning from the CS state is obtained numerically. The discontinuous 
increase of the depinning threshold at the critical point is a peculiarity of the piecewise pinning potential 
and is replaced by a sharp, but continuous rise for other pinning potentials. 



from the coherent phase at weak disorder to the incoherent phase at strong disorder is discontinuous 
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FIG. 3: Coherence (right) and velocity (left) as functions of drive for ho/K = 1 (red curves, showing the 
CS — * CM continuous depinning) and ho/K = 4 (black curves, showing the IS — > CM hysteretic depinning). 
The curves are obtained by numerical integration of Eq. (|15|1 as F is ramped up and then down. 

for piecewise parabolic pinning. It also becomes hysteretic for most other pinning potentials. [~LtJ| 
When the driving force is ramped up adiabatically from the coherent state, the system depins 
continuously at F* to a unique sliding state and the coherence grows smoothly towards unity. 
When the force is ramped up from the incoherent static state, where r = 0, the system depins 
discontinuously at F| = F sp , with F sp = ho the single particle threshold, given by the maximum 
pinning force. This is easily understood from Eq. H15|) as when r = the displacements can remain 
decoupled even for finite values of K. At Fl the coherence also jumps discontinuously to a finite 
value, as shown in Fig. |3J When the force is ramped back down from the sliding state, the system 
repins at the lower force F\, where the coherence also jumps back to zero. The moving state 
is always coherent in mean field, although incoherent sliding states should be possible in finite 
dimensions [l^ |. The origin of the mean-field hysteresis observed when the system depins from a 
static incoherent state is easily understood. When the force is ramped up each degree of freedom 
depins essentially independently. Once the medium starts sliding, disorder becomes less important 
and the system becomes coherent and therefore much stiffer. Such a system cannot easily adjust 
to disorder and therefore when the force is decreased remains in a sliding state down to a lower 
force. This type of hysteresis is clearly strongly enhanced in mean-field, where only states with 
zero or perfect coherence are possible. Hysteresis has been observed in numerical simulations of 
the PS model in three dimensions [3], although this work did not establish conclusively that the 
hysteresis survives in the limit of infinite system size. Preliminary work by us [3| suggests that 
there is no hysteresis in infinite systems in three dimensions. 
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IV. CONCLUSION 

We have discussed the depinning transition of two classes of models that allow for history- 
dependent response. In mean field theory both models exhibit a tricritical point as a function 
of disorder strength. At weak disorder depinning is continuous and the sliding state is unique. 
Above the tricritical point depinning is discontinuous and hysteretic. Numerical studies of these 
models are currently underway to establish how much of the mean field behavior survives in finite 
dimensions. Simulations of the PS model in three dimensions suggest that the hysteresis disappears 
for large system sizes In this model the hysteresis may indeed be an artifact of 

the mean field approximation that yields a sharp distinction between coherent and incoherent 
static states. In finite dimension this may be replaced by a smooth growth of local coherence 
over a broad range of length scales, without discontinuous jumps. The mean-field tricritical point 
may then become a strong crossover between continuous depinning transitions characterized by 
(nonuniversal) exponents (3 < 1 at weak disorder to j3 > 1 at strong disorder j^j. Hysteresis seems 
more robust in the VE model, or in general in models with local inertia. Preliminary numerical 
simulations of an anisotropic version of the VE models in two dimensions indicate a stubborn 
persistence of hysteresis with increasing system size even for relatively weak disorder. In this case 
the existence of hysteresis may not be a good test of the precise nature of the depinning transition 

Q. 

Once local inertia or topological defects are introduced in the model, various depinning scenarios 
are possible. The depinning may be discontinuous with hysteresis (like an equilibrium first order 
phase transition) or with hysteresis that vanishes in the infinite system limit. Another possibility 
is that the transition in continuous and critical, in the sense that it is possible to identify diverging 
correlation lengths as the depinning threshold is approached adiabatically from above or from 



291 ] and more recently explored 



below. Finally, as originally suggested by Ramanathan and Fisher 
by Maimon and Schwarz 2£|, even critical behavior with hysteresis that survives in the infinite 
system limit is possible. Sorting out this various scenarios for the models discussed will require 
more extensive numerical studies in finite dimensions. 
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